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INTRODUCTION 


Fast  Fourier  Transform  (FFT) 


An  external  subroutine  is  available  which  performs  the  discrete 
Fourier  Transform  on  a multidimensional  array  of  floating  point  data. 
The  data  may  be  either  real  or  complex.  The  transform  values  are 
always  complex  and  are  returned  to  the  array  used  to  carry  the  original 
data.  The  dimensionality  of  the  input  array  is  limited  only  by  machine 
storage  capacity.  The  number  of  input  samples  is  not  limited  to  a 
power  of  two  and,  in  fact,  can  be  prime  in  any  or  all  dimensions.  The 
subroutine  is  called  by  a one  line,  six  argument  call  statement  in 
the  user  program. 

v 

Procedure  for  FFT 

The  FFT  subroutine  is  called  with  the  following  statement: 

CALL  FOURT(DATA,NN,NDIM,ISIGN,IFCRM,WQRK) 


♦The  Program  described  in  this  memo  has  been  adapted  for  the  UNIVAC  1108 
by  the  authors  from  a Lincoln  Laboratory  Technical  Note  "Three  Programs 
That  Perform  the  Cooley-Tukey  Fourier  Trams form"  by  N.  M.  Brenner  on 
28  July  1967.  The  AD  No.  is  AD  657019. 
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where  j 

UATA  « An  N-dlmenslonal  complex  array  used  to  told  the  real  and 
imaginary  parts  of  the  Input  data,  and  the  real  and 
imaginary  parts  of  the  transformed  values . The  real  and 
imaginary  parts  of  a datum  must  be  stored  In  adjacent  core 
locations  as  done  by  the  complex  array  declarer  in  fortran  V. 
For  real  input  data,  the  array  DATA  is  a complex  array  with 
the  Imaginary  part  set  to  0.0. 

NN  * A single  dimensional  integer  array  equal  in  length  to  the 

number  of  dimensions  in  the  transform.  The  first  element  in 
the  array  is  equal  to  the  length  of  the  first  data  dimension, 
the  second  element  is  equal  to  the  length  of  the  second  data 
dimension,  etc. 

NDXM  ■ The  number  of  dimensions  used  for  the  array  DATA. 
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1FGRM 


1-1  for  a 
1 + 1 for  ai 


a forward  transform, 
an  Inverse  transform. 


+ 1 if  data  and  transform  are  complex. 

0 is  data  is  real  (zero  Imaginary  values). 

If  it  is  "0"  the  Imaginary  parts  of  the  data  should  be  set 
to  zero.  (Transform  will  be  complex). 

"0"  if  all  dimensions  are  equal  to  a power  of  two.  If  not  it 
must  be  a single  dimension  floating  point  array  equal  to  twice 
the  length  of  the  largest  data  array  dimension  not  a power  of 
two. 


Special  Notes; 

I.  There  are  no  error  messages,  error  halts  or  error  returns  in  this 
subroutine.  If  NDD4  or  any  Nn'('i)  is  less  than  cxie,  the  program  returns 
immediately. 

?.  The  data  is  assumed  to  form  one  cycle  of  a periodic  function.' 

J.  The  data  is  assumed  to  be  equispaced  in  each  dimension1?  If  the 
equispaced  interval  Is  PT,  the  resulting  transform  will  be  equispaced 
from  0 to  2 Tl  (M-l)/(NDT).  The  upper  limit  la  identified  with  II /(NIT ) 
and  all  points  above  the  foldover  frequency  11  /NOT  are  identified  with 
the  corresponding  negative  frequency. 
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4.  If  an  Inverse  transform  is  performed  on  an  array  of  transformed 
data,  the  original  data  will  reappear  multiplied  by  N1*N2*  . . . 

5.  The  running  time  is  extremely  fast  if  the  dimensions  are  rich 
in  factors  of  two.  The  particular  running  time  is  dependent  on  the 
factorization  of  each  dimension.  Approximate  times  for  a one 
dimensional  complex  transform  are  as  follows: 

N Factorization  Time  for  Complex  Transform  (SEC) 


N Factorization  Time  for  Complex  Transform  (SEC) 

4094  2 x 23  x 89  **0 

4096  32  x 5 x 7 x 13  12 

4096  2 12  3.1 

4097  17  x 241  90 

4098  2 x 3 x 683  2**0 

4099  PRIME  1430 

4100  22  x 52  x 4l  20 

6.  For  the  use  of  real  data  input  the  designation  IF  CRM  ■ 0 can 
reduce  running  time  by  as  much  as  forty  percent. 

7.  A user  must  allow  about  1750  core  locations  for  the  instructions. 

8.  The  source  and  relocatable  elements  are  on  File  2 of  Cur  Tape  U-183. 
Examples  of  FFT  Program 

A.  One  dimensional  forward  transform  of  real,  data  of  length  64. 


COMPLEX  D (64) 
Do  10  I = 1,  64 
10  D(I)  * (R,  0.0) 


6 Declare  Data  Array 
@ Zero  Imaginary  Values 


Fortran  Program 


CALL  FOURT  (D,  64,  1,  -1,  0,  0) 
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Data  Array- 
Length  of  Array 
One-Dimension 
Forward  Transform 
Data  is  real 
No  working  array  needed 


B.  Three-dimensional  forward  Fourier  transform  of  a complex  array 
dimensioned  32  by  25  by  13. 

DIMENSION  WORK  (50),  NN  (3) 

COMPLEX  D (32,  25,  13)  @ Declares  Data  Array 

DATA  NN/32,  25,  13/  @ Supplies  Dimension  Info  to  NN 


CALL  FOURT  (D,  NN,  3,  -1,  1,  WCRK) 


Data  Array 

Integer  Array  of  Length  3 to  Supply  Dimension  Info 
Three  Dimensional  Transform 
Forward  Transform 
Complex  Transform 

Array  of  size  50  (2  x 25)  since  25  is  largest  non  ] 
of  two  dimension. 


WCRK 


Conclusions 


The  discrete  Fourier  Transform  presented  is  extremely  flexible 
in  nature.  Because  of  this  flexibility  considerable  storage  is 
required  compared  to  some  special  purpose  FFT's.  It  is  felt  that 
the  main  use  of  this  FFT  should  be  as  cm  analytical  tool  as  opposed 
to  a "real  time''  processing  technique. 
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